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Решена задача о промерзании грунта с использованием математической модели, существенно отличающейся от классической 
постановки Стефана, без явного выделения границы фазового перехода. Задача решена методом конечных разностей с приме- 
нением неявной итерационной разностной схемы. Установлено, что предложенный подход приводит к результатам, отличаю- 
щимся не более чем на 4,4 % по температурам и по координате фронта промерзания при решении одномерной задачи. 


Многие процессы теплообмена связаны с изме- 
нением агрегатного состояния материала. При этом 
теплофизические характеристики материала изме- 
няются скачкообразно, и при фазовых переходах 
происходит выделение или поглощение теплоты [1]. 
Решение подобного рода задач имеет большое прак- 
тическое значение в металлургии, строительной те- 
плотехнике и в других прикладных дисциплинах. 

В большинстве случаев решения задач о движе- 
нии границы между жидкостью и твердым телом 
проводят по методу Стефана [2]. Следует отметить, 
что данный метод применим для одномерной по- 
становки задачи, но значительно сложнее обстоит 
дело с двухмерной или с трехмерной задачей. Это 
является недостатком метода Стефана. В данной 
работе рассматривается другой подход к моделиро- 
ванию фазового перехода, отличный от классиче- 
ской задачи Стефана [2]. 

Рассмотрим этот подход на примере задачи о 
промерзании насыщенного влагой грунта. Анализ 
проводится на примере пластины бесконечной дли- 
ны, толщиной Ь с некоторой начальной фиксиро- 
ванной температурой Г 0 . В начальный момент вре- 
мени на поверхности грунта устанавливается тем- 
пература Т„, ниже температуры замерзания Г 3 воды. 
В результате грунт начинает промерзать. Образуется 
промерзший слой переменной толщины <*, завися- 
щей от времени. Подвижная граница промерзания 
всегда имеет температуру кристаллизации. На этой 
границе происходит переход из одного агрегатного 
состояния в другое, при этом выделяется теплота 0. 

Для описания процесса охлаждения грунта при 
промерзании используется модель [3], в которой яв- 
но не выделяется граница промерзания в граничных 


условиях и область решения не разбивается на две 
подобласти, соответствующие разным агрегатным 
состояниям. Процесс тепловыделения при кристал- 
лизации описывается Источниковым членом в ура- 
внении теплопроводности. Значение объемной ско- 
рости промерзания \Ѵ рассчитывается по скорости 
продвижения изотермы Т=Т Г Теплофизические 
свойства материала считаются постоянными. Задача 
рассматривается в неподвижной системе координат. 


Математическая постановка задачи в рамках 
рассматриваемой модели имеет вид: 


ВТ д 2 Т О ■ \Ѵ 

т > 0, 0< х< Ь, 


( 1 ) 


где с - теплоемкость; р - плотность; Т- температу- 
ра; Я - коэффициент теплопроводности; г - время; 
х — координата; I - толщина слоя грунта; (2 3 - тепло- 
та кристаллизации; \Ѵ=р ■ V— массовая скорость про- 
мерзания; V- линейная скорость кристаллизации. 


Начальные условия: 

г=0, Т=Т 0 , 0 <х<Ь, (2) 

т=0, <^=0. (3) 

Т 0 - начальная температура. 

Граничные условия: 

х=0, Т=Т„, (4) 

Х=%, Т=Т= СОП8І, (5) 


х = Ь, 



( 6 ) 


Нелинейное дифференциальное уравнение те- 
плопроводности для слоя грунта с граничными 
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условиями первого рода на поверхности охлажде- 
ния решена методом конечных разностей [4] . Раз- 
ностные аналоги исходных дифференциальных 
уравнений решены методом итераций с примене- 
нием неявной четырехточечной разностной схемы. 
Область решения задачи массопереноса имела по- 
движную границу, поэтому вычисления проводи- 
лись на нерегулярной сетке, учитывающей продви- 
жение фронта промерзания вглубь грунта с пере- 
менной во времени скоростью. Особенность пред- 
ложенного метода основывается на том, что ско- 
рость промерзания рассчитывается по температу- 
рам в узлах разностной сетки. С течением времени 
толщина промерзания увеличивается. Численное 
значение Е, рассчитывается по формуле: 

= + , (7) 


здесь к - шаг разностной сетки по х; і - номер узла 
сетки, в котором значение температуры достигло 
температуры замерзания воды; 7]_ ь Т, - значение 
температуры в данный момент времени в (/-1) и (і) 
узлах сетки. 

Линейную скорость кристаллизации можно вы- 
числить, используя следующее выражение: 


ѵ _ %(т+Ат)-%(т) н 
Аг 


( 8 ) 


где ^(г+Дг) и <^(г) — толщина промерзшего слоя в 
рассматриваемый и предыдущий моменты време- 
ни; Аг - шаг разностной сетки по времени. 

Поскольку поставленная задача решалась в од- 
номерной постановке, то данный метод был проте- 
стирован по задаче Стефана о промерзании насы- 
щенного влагой грунта. Решение Стефановской за- 
дачи проводилось численно, по выше указанному 
алгоритму (1-8). Было выделено две подобласти 
решения. На каждом шаге по времени определя- 
лась толщина промерзания (граница раздела обла- 
стей), а затем рассчитывалось температурное поле. 

Температура замерзания различных грунтов зави- 
сит от влажности и консистенции, для водонасыщен- 
ных песчаных, текучих и текуче-пластичных глини- 
стых грунтов данная температура близка к 273 К [5]. В 
этом случае задачу о промерзания грунта можно све- 
сти к задаче образования льда в воде [2]. Поэтому при 
проведении расчетов использовались значения: 
С=0,03 м; 7’ 0 =247 К; р = 920 кг/м 3 ; р=1000 кг/м 3 ; 
с, =2 100 ДжДкг-К); с=4200 Дж/(кг-К); А=2,7 Вт/(м-К); 
А, =0,68 ВтДм-К); 7=273 К; 0=334 кДж/кг [6, 7]. 

На рис. 1 приведено типичное распределение 
температур по х при условии Т=212 К. 


На рис. 1 представлены зависимости температу- 
ры от координаты, полученные при решении зада- 
чи о промерзании грунта двумя различными подхо- 
дами. Кривая 1 соответствует выше предложенной 
модели (1-8), а кривая 2 - модели Стефана. Как 
видно, температурные поля для разных решений 
незначительно отличаются друг от друга. На рис. 2 
показано изменение величины промерзшего слоя 
от времени. При расчете толщины промерзшего 
слоя Е, на момент времени 200 с, были получены 
следующие результаты: 1,34 мм для решения пер- 
вым способом, 1,28 - вторым, рис. 2. Анализ полу- 
ченных результатов показывает, что относительная 
погрешность вычислений величины Е, составля- 
ет 4,4 %. Это позволяет сделать вывод о целесооб- 
разности использования модели (1-8) при реше- 
нии задачи с учетом теплоты кристаллизации. 
Примером задач с фазовым переходом также явля- 
ются такие процессы как конденсация и плавле- 
ние. Следовательно, выше предложенная модель 
может быть применима и для их решения. Этот 
подход позволяет упростить процедуру расчета ма- 
тематического моделирования многомерной зада- 
чи с фазовым переходом. Следует отметить, что для 
проведенного исследования был использован от- 
носительно небольшой перепад температур. 


т,к 



Рис. 1. Распределение температуры по координате х: 1) мо- 
дель (1-8), 2) задача Стефана [2] 
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Рис. 2. Зависимость толщины промерзания от времени: 
1) модель (1~8), 2) задача Стефана [2] 
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